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O-t- Abstract 

i We present a new method allowing us to derive the long-time and large-distance asymp- 

. totic behavior of the correlations functions of quantum integrable models from their exact 

representations. Starting from the form factor expansion of the correlation functions in 
finite volume, we explain how to reduce the complexity of the computation in the so-called 
interacting integrable models to the one appearing in free fermion equivalent models. We 
^ . apply our method to the time-dependent zero-temperature current-current correlation func- 

' tion in the non-linear Schrodinger model and compute the first few terms in its asymptotic 

. expansion. Our result goes beyond the conformal field theory based predictions: in the time- 

QP ' dependent case, other types of excitations than the ones on the Fermi surface contribute to 

^—^ the leading orders of the asymptotics. 

1 Introduction 

The successful resolution of many-body quantum integrable models through Bethe Ansatz [31 
^ ■ EH |83l EH] naturally raised the question of the computation of their correlation functions. In 

^ I its full generality, this problem appeared first as extremely complex due to the complicated 

combinatorial representation of the eigenfunctions. However, for the particular values of the 
coupling constants for which the models are equivalent to free fermionic ones, the use of specific 
methods (such as the Jordan- Wigner transformation) led to numerous explicit results [661 16^ 
IMlESEQliaiTellHilHIHSSQllTalElllSa in this context, different 

types of effective representations have been obtained for the correlation functions, notably in 
terms of determinants. It has been possible in particular to express two-point functions in 
terms of Fredholm determinants of compact operators I + V, where the integral kernel V 
belongs to the class of integrable integral operators |64| [65| \57\ \T7\ 116) . These determinants 
representations have been used to compute the asymptotic behavior of the correlations functions 
at large distances [301 [HI EHl ESI [Ml H • 
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It was the understanding of the underlying algebraic structures in truly interacting integrable 
models which finally opened the way to the computation of their correlation functions. In 
particular, the algebraic version of the Bethe Ansatz (ABA) |24) enabled Izergin, Korepin and 
their collaborators in Leningrad to obtain the first exact representations for the correlation 
functions of the XXZ and the non-linear Schrodinger models away from their free fermion 
points (see |56| and references therein). These representations, based on finite-size determinant 
representations for the norms of Bethe eigenstates |25| [26l [55] and for more general scalar 
products |82| I56| . were however not completely explicit due to the introduction of dual fields 
to handle the combinatorial difficulties of the ABA approach. The first explicit representations 
(multiple integral representations for elementary building blocks of the XXZ chain at zero 
temperature) were obtained in the 90' by Jimbo, Miki, Miwa and Nakayashiki through q- 
vertex operators and solutions of the q-KZ equation |37| I^Uj I38| . A few years later, these 
representations, together with their generalization to the case of an external magnetic field, were 
reproduced by the Lyon group in the ABA framework thanks to the resolution of the so-called 
quantum inverse scattering problem |53| [5l] . The Lyon group was also able to provide finite-size 
determinant representations for the form factors of the finite XXZ chain |53| . Developments of 
the latter approach led to integral representations for the two-point function |50[ l52| , to their 
generalizations in the temperature case |27| |28] , in the time-dependent case |51 ] , and in all these 
cases taken together [80j- On the other hand, new interesting algebraic representations for the 
correlation functions were also obtained recently using reduced q-KZ equations [9l [TOl H] , which 
led to the discovery of a hidden Grassmann structure in the XXZ model |12^ |8l [131 [^T| [Sl [7] . 

In this framework, one of the most challenging problem is the derivation, from these ex- 
act representations, of the long-distance asymptotic behavior of the correlation functions. In 
general, zero-temperature correlation functions of models with a gapless spectrum are expected 
to decay algebraically at large distances. The exponents governing this power-law behavior, 
called critical exponents, are believed to be universal, i.e. not depending on the microscopic 
details of the model, but only on its overall symmetries. Predictions concerning these exponents 
were obtained via the Luttinger liquid approach |69| [29} [30] I31| or the connection to Conformal 
Field Theory |79[ [Hj [H [151 E| , together with exact computations of finite size corrections to the 
spectrum |20l I86| [2T 1 187 | [22]. However, until very recently, it was not possible to confirm these 
(indirect) predictions via some direct computations based on exact results for an integrable 
model out of its free fermion point. 

To tackle this problem, one may notice that, in the ABA approach, the correlation func- 
tions can be written in terms of series of multiple integrals that can be though of as some 
multidimensional generalizations of the Fredholm determinants which appear in the study of 
the free fermionic models. One may therefore attempt to extend and adapt the large arsenal of 
techniques available for free fermionic correlators to models away from their free fermion points. 
The main obstacle here comes from the highly coupled nature of the various series of multiple 
integrals representations for the correlation functions. Actually, until two years ago, there were 
only few cases out of the free fermion point where the analysis could have been performed to 
some extend [48 [ I49 | [6T] . A major breakthrough was achieved in |45) . which proposed a system- 
atic method to perform the asymptotic analysis of the long-distance asymptotic behavior of the 
correlation functions of densities in integrable models described by a six-vertex type i?-matrix 
and solvable by the ABA. There, the idea was to extract the asymptotic behavior of the series of 
multiple integrals representing the correlator by making some connections with the asymptotes 
of Fredholm determinants. Although fully functional, the method involved several subtle steps, 
and in particular a quite complex summation process. The final result itself, nevertheless, hap- 
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pened to be quite simply given in terms of a properly normalized form factor |45t 146) . This fact 
naturally suggests that the method proposed in [45] might be simplified if one works directly 
at the level of the form factor expansion of the correlation functions. 

It is the aim of the present article to expound such a simpler and more direct approach, 
based on the form factor expansion of the correlation functions in the ABA framework, on the 
finite-size determinant representation for the form factors in finite volume and on some of their 
main features in the thermodynamic limit. The main advantage compared to |45| is that it 
relies on the study of more natural physical objects: whereas in pS] the correlation function 
was written as a series over some elementary objects expressed in terms of bare quantities (such 
as energy, momentum), the so-called cycle integrals, here we chose to pursue our study of form 
factors a step further, and to use some of their thermodynamic properties (and in particular their 
representation in terms of dressed physical quantities) from the very beginning. This results into 
a major simplification of the summation process. Indeed, in |45| . the asymptotic summation 
of the series of cycle integrals was quite subtle. In particular, it involved a summation of a 
whole class of corrections, sub- leading to each term of the series, but nevertheless contributing 
to the leading term in the final result: they happened to be responsible, once upon summation, 
for the dressing of bare physical quantities. Since in our present approach the series over form 
factors is already written in terms of the dressed quantities, the summation is much simpler: 
the series can be connected to a Fredholm determinant and the leading asymptotic behavior 
of the correlation function is directly given by the leading asymptotic behavior of this Fredholm 
determinant. No subtle identification of sub- leading terms is needed for the result. 

Let us be slightly more precise about the difi^erent steps of our method. Its starting point is 
the form factor expansion of the correlation function in finite volume. Using the determinant 
representation of the form factors in finite volume |8HI53j. one can analyze their thermodynamic 
limit in the particle/hole picture \82\ H6l H3] . and decompose their finite-size expression into 
some (universal) "singular" or "discrete" part, which is essentially free-fermion-like and contains 
the whole non-trivial scaling behavior of the form factor in the thermodynamic limit, and some 
(model-dependent) "regular" part, which admits a smooth thermodynamic limit |43) . The idea 
is to reduce the summation over form factors to a summation over their finite-size discrete parts, 
the (thermodynamic limit of the) regular parts being treated as some dressing functionals. The 
main problem is that, for interacting models, the series we consider is highly coupled. The 
origin of this coupling is twofold: on the one hand, the regular part is a complicated function 
of the particle/hole Bethe roots; on the other hand, the discrete part can be seen as a certain 
functional of the shift function between the considered excited state and the ground state, and 
hence introduces non-trivial couplings between the particle and hole-type Bethe roots. It is 
possible to slightly simplify the structure of the different terms by neglecting, on the level of 
the finite-size series, some contributions that should vanish in the thermodynamic limit. We 
thus obtain an effective (finite-size) form factor series with a simpler particle/hole dependence 
than in the original series. This series can be linked to a decoupled one, reduced to the sum of 
the discrete parts of the form factors with a shift function that does not depend on the excited 
state, exactly as in a (generalized) free fermionic model. The sum can then be computed and 
recast in terms of a determinant which tends to a Fredholm determinant in the thermodynamic 
limit. The leading asymptotic behavior of the correlation function then follows directly from 
the leading asymptotic behavior of the Fredholm determinant |60| . 

We stress that, within this method, the computation of the asymptotic expansion of the 
correlation functions for "truly interacting" models and for their free fermionic limits carry 
almost the same complexity, in as much as we offer a way of understanding the interacting case 
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as a smooth deformation of the free fermionic one. This joins in some sense the spirit of the 
works [12J IHl [T3t Wi\ [6l [7] on the hidden Grassmann structure in the XXZ model with the main 
difference that, in our method, the free fermionic structure is only used as an intermediate step 
so as to carry out certain algebraic manipulations in a simpler way and deform the range of 
summations to a kind of analog of the steepest descent contour. 

We choose to expose our method in the framework of the quantum non-linear Schrodinger 
model, which is probably the best understood and simplest integrable interacting model. Note 
however that, with minor modifications specific to the structure of the model under investiga- 
tion (more complex structure of the configuration of the pertinent Bethe roots above the ground 
state, the appearance of trigonometric functions . . . ), this method is in principle applicable to 
a very wide class of integrable models whose form factors are known and admit a finite-size 
determinant representation such as those obtained in |8 H [77]. We apply our method to the 
derivation of the first few terms in the long-distance and large-time asymptotic expansion of 
the correlation function of currents at zero temperature (see Section [3]). The case of the re- 
duced density matrix, together with some more precise mathematical details about the method, 
will be considered in a further publication |59) . We will see that the asymptotes in the time- 
dependent case are not only issuing from excitations in the vicinities of the Fermi boundaries 
(which correspond to the region of the spectrum taken into account by the CFT/Luttinger 
liquid approach), but also from the excitations around the saddle-point Aq of the "plane-wave" 
combination xp{X) — te{X), where p and e are respectively the dressed momentum and energy of 
the excitations. More precisely, the exponents governing the power-law decay of the correlation 
function are given in terms of the different shift functions between the ground state and excited 
states having one particle and one hole either on the Fermi boundaries ibg or at the saddle- 
point Aq- The associated amplitudes are given by the corresponding, properly normalized, form 
factors of currents. 

The paper is organized as follows. In Section [21 we introduce the non-linear Schrodinger 
(NLS) model and the different notations, relative to the description of the space of states and 
the thermodynamic limit of the model, that will be used throughout the article. This enables 
us, in Section[3l to present our main result: the leading large-distance and long-time asymptotic 
behavior of the (zero-temperature) correlation function of currents. The remaining part of the 
article is devoted to the derivation of this result. We establish in Section!?] a (finite-size) form 
factor series representation for the correlation function (or more precisely for its generating 
function). In Section [5l we explain how to relate our highly coupled series to a decoupled, free- 
fermion type series that can be summed up into a finite-size determinant. Finally, in Sectional 
we take the thermodynamic limit of the previous result, hence representing the form factor 
series in terms of a Fredholm determinant, and explain how to derive the first leading terms 
of the long-time/large-distance asymptotic behavior of our series from those of the Fredholm 
determinant. Technical details are gathered in a set of Appendices: explicit representations for 
form factors are given in Appendix lAj the summation of the free-fermion type series is performed 
in Appendix [Bl the notion of functional translation operator, which is used as a technical tool 
to link our series to a decoupled one, is introduced in Appendix [Cl in Appendix [D1 we show 
how to relate the summed-up form factor series to the type of series that was obtained in [45] by 
expanding the master equation; finally, in Appendix |E1 we explain how to control sub-leading 
corrections in the asymptotic series using the so-called Natte series introduced in |60| . 
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2 The NLS model: notations and definitions 



The quantum non-hnear Schrodinger model is described in terms of quantum Bose fields ^{x,t) 
and obeying to canonical equal-time commutation relations, with Hamiltonian 

L 

H = j ^d^^^d:,<If + c^^^^<If^-h^''^^dx. (2.1) 



The model is defined on a finite interval of length L (we take the thermodynamic limit L — t- oo 
later on) and is subject to periodic boundary conditions. In (j2.ip . c denotes the coupling 
constant and h the chemical potential. In this article, we focus on the repulsive regime c > 
in presence of a positive chemical potential h > 0. 

Due to the conservation of the number of particles, this quantum field theory model is 
equivalent to a (quantum mechanical) one-dimensional many-body gas of bosons subject to 
point-like interactions. Such a model was first introduced and solved by Lieb and Liniger in 
1963 [681 I67| . where it was used as a test for Bogoliubov's theory. It was then thoroughly 
investigated. In particular, in 1982, Izergin and Korepin |36] (see also |56| ) proposed a lattice 
regularization of the model allowing to implement the ABA scheme [24] in finite volume. 

In this article, we consider the zero-temperature time-dependent correlation function of cur- 
rents, i.e. the ground state expectation value ( j(x, i) j(0, 0) ), where j{x,t) = ^^{x,t)"^{x,t). 
Our aim is to evaluate, in the thermodynamic limit L — )• oo of the model, its large-distance and 
long-time asymptotic behavior. Prior to writing down our result, we introduce some necessary 
notations. Namely, we describe the space of states of the model and provide definitions of 
several thermodynamic quantities which appear in our result. 

In the algebraic Bethe ansatz framework, the eigenstates | ■ip{{X}) ) of the Hamiltonian (j2.ip 
are parametrized by a set of spectral parameters Ai, . . . , Aat solution to the system of Bethe 
equations. In their logarithmic form, these Bethe equations read 

N 

Lpo{\j) + Y,^{\j-\k) = 2TTnj, j = l,...,iV, (2.2) 

k=l 

in terms of a set of integer (for odd) or half-integer (for A^ even) numbers ni, . . . ,njv. The 
bare momentum po{\) and the bare phase 6{\) appearing above are given by 

Po(A) = A, 0(A) = i log (^). (2.3) 

\ic — A/ 

The energy of the state | V'({'^}) ) is 
TV 

^({A}) = J]eo(A,), with eo{X) = X^-h. (2.4) 

i=i 

The solutions to the system of Bethe equations have been studied by Yang and Yang |89| . 
In particular, it has been proven that all solutions are sets of real numbers and that a set of 
solutions {Aj} can be uniquely parametrized by a set of (half-)integers {n^}. 

In the following, Aj, j = 1, . . . , A^, will denote the Bethe roots describing the ground state of 
the Hamiltonian (12. Ih . They correspond to the solution to (12. 2p associated with the following 
choice of (half)-integers: Uj = j — {N + l)/2, where the number of particles A^ is fixed by the 
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value of the chemical potential h. In the thermodynamic limit (L — )• oo, N ^ oo with N/L 
tending to some finite average density D), the parameters Xj densely fill a symmetric interval 
[—q, q] of the real axis (the Fermi zone) with some density distribution p(A) solving the following 
linear integral equation 

pW -^Jm-f,) p(^) d^ = (2.5) 
where K denotes the Lieb kernel 



A2 + C 



Excited states of (12. ip correspond to other solutions of the Bethe equations. For technical 
purpose, it is actually convenient to consider the twisted Bethe states \ ^k({/^}) ) instead. These 
are parameterized by solutions /x^^, . . . of the twisted Bethe equations, 

iPo(/U^,) + X;^(^£, -/^4) = 2vr(£, -^^^)+i/3, j = l,...,N^. (2.7) 
k=l 

Here (3 is some purely imaginary parameter, k = e^, and < ^2 < " " " < ^A^^ are some integers. 
Note that, the set of roots of (12. 7p being completely defined by the set of integers ij, we have 
chosen to label them accordingly. Eigenstates of the Hamiltonian (|2.ip are obtained as limits 
when K — )• 1 of the twisted Bethe states | ^^^{{fi}) ). 

For purely imaginary f3, the roots of (j2.7p are real. The state parameterized by the solutions 
of (12. 7p with ij = j , j = 1, . . . ,Ni^, is called the K-twisted ground state in the A^^-sector. In 
the thermodynamic limit, excitations above this k- twisted ground state correspond to solutions 
such that most of the ij^s coincide with their value for the ground state: ij = j except for 
n integers (with n remaining finite in the N,^ — )• 00 limit) hi, . . . ,hn G {!,..., A^^^;} for which 
ih^ = Pk ^ {1, • • • , Nk}- The integers hk represent "holes" with respect to the distributions of 
integers for the ground state in the A^^-sector, whereas the integers pk represent "particles". 
For such an excited state | V'K({^£j}))j the counting function 

Uc) ^ Uu I {/.,J) = ^poiu) + ^ 0(^ - + - (2.8) 

k=l 

which is monotonously increasing, defines unambiguously a set of real "background" parameters 
lij, j G Z, as the unique solutions to CK(/^i) = j/^- Among this set of parameters, the N,^ 
solutions of the twisted Bethe equations with integers ij correspond to the solutions {fJ-lj ) = 
ij/L, j = 1, . . . ,Ni^, and the particle rapidities fip^ (respectively the hole rapidities ^/i^), a = 
1, . . . ,n, correspond to the solutions ^^(/^pa) — Pa/L (respectively inil^ha) — ^a/L). 

Due to the particular type of correlation function we consider in this article (current-current 
correlation function) , we will restrict our study to states having the same number of parameters 
as the ground state, i.e. Nf^ = N. To describe the position of the roots of an excited state 
I i^Kiiptj}) ) with respect to the ground state roots {Aj}, it is convenient to define the (finite- 
size) shift function 

F{u) = F(a;|{^,J) = L[^{lo) - U^)] = [^('^ " ^'^) " ^('^ " ^4-)] + ^> (2-9) 

k=l 
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in which is the excited state's counting function (|2.8|) whereas ^ denotes the ground state 
counting function (at k = 1). The spacing between the root \j for the ground state and the 
background parameters define^l by S,K{lJ-j) = j /L is then given by 

- = + 0{L-^), j = 1, . . . , A^, (2.10) 

in which F{\) is the thermodynamic hmit of the shift function, solution of the integral equation 

^ . n 

FW - ^ / i^(A - /i) F{^^) d/z = ^ - [^(A - ^,,,) - e{\ - ^,J] . (2.11) 

4 k=i 

Introducing the dressed phase 0(A, /i) and the dressed charge Z{\) as the respective solutions 
of the linear integral equations 

•/-(A,//)-^ / ^(A-^)'A(c^,/^)du; = ^^^^, (2.12) 
-q 
q 

Z{X) - ^ y i^(A - /x) Z(^) dfs = 1, (2.13) 

-a 

we get 

■o ri 

F{\) = '£Z{X) - [<^(A,/ipJ - HK^^hJ]■ (2.14) 

a=l 

Other important thermodynamic quantities that we need to introduce in order to formulate 
our result are the dressed energy e(A) and the dressed momentum p(A), defined as 

q 

e{X) - — I K{X- fi)e{fi)dfi = eo{X), with e{±q) = 0, (2.15) 
27r J 

-q 

q A 
p(A) = po{X) + ( e{X- is) pill) dii = 2n [ pill) dp. (2.16) 



-g 

Note that expression (|2.16p enables us to express, in the large L limit, the counting function 
(12. 8p of any excited state with a finite number of particles and holes above the A^-particle ground 
state in the following form: 

U^\{p,^})=Cicj) + OiL-'), with ^(o;) = i-p(L.) + |, (2.17) 

D being the average density in the thermodynamic limit iN/L D). In particular this means 
that, in the thermodynamic limit and up to 0(1/L) corrections, the counting function does not 
depend on the particular localization of the corresponding Bethe roots {pgj}: it is the same for 
all particle/hole- type excited states. 



^Since the counting function depends on the particular excited state we consider, so does the set of background 
parameters: these are therefore defined for each excited state {fie } separately. 
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3 The main result: large-distance/long-time asymptotic behav- 
ior of the correlation function of currents 



We are now in position to formulate our main result, namely the leading asymptotic behavior of 
the zero-temperature correlation function of currents j(0, 0) ) in the large-distance and 

long-time regime, with x/t = const. 

Let u{X) = p{X) — ^£{X), where p is the dressed momentum (I2.16P and e the dressed 
energy (|2.15p . the ratio t/x being fixed (and non-zero). We assume that this function has a 
unique saddle-poini0 Aq on R (i.e. u' admits a unique zero Aq on M, which moreover satisfies 
u"{\q) < 0). One distinguishes two different regimeq^ according to whether Aq G] — q, q[ (time- 
like regime) or Aq ^ [— q\ (space-like regime). Then, at large-distance and long-time (x — )• +oo, 
t — )• -|-oo, x/t = const) and in the space-like regime, 



(,(x,t),(0,0)) = (^) a + 

2cos(22;pir) • \T'^ ' 



-<?i 



-i{x — vpt)]^'' [i{x + vpt)] 



Z2 



(l + o(l)) 



^ Vte"(Ao) - xp"iXo) [-i(x - t;^t)][^.'°W-i]' + ^^^^^[F,">i-q)P + ° 
, ^/2^e-ip-(Ao) ei^\Pi^o)+PF]-ueiXo) . \:pXo^f ^ 

whereas in the time-like regime, 

(,(x,t),(o,o)) = (f) -^jj^^i^ + oiD) 



2cos{2xpf) ■ \ j^qf 
-i{x — vpt)]^^ [i{x + vpt)]^'^ 



(l + o(l)) 



^ yt?W^W^H(x-t;^t)][^'o(^)+i]^ [^{x + vpt)\n^^-'^^^' ^ 

^ v*^"(Ao) - xp"(xr) [_i(x - z;^t)][^A7wi^ [i(^ + v^t)f^:(-''^-'y' ^ ^^'^^ 

In these expressions, pi;' = zizp{zizq) = ttD is the Fermi momentum, vp = £'iq)/p'{q) is the 
Fermi velocity, whereas Z = Z(ibg) is the value of the dressed charge on the Fermi surface. 
The exponents of the last two terms in (I3.ip - (l3.2p are given in terms of the values at the Fermi 
boundaries itg of some specific shift functions (|2.14p Fy^^ (at /3 = 0) between the ground state 
and an excited state with one particle at fj,p and one hole at Hh- the shift functions 

Fi«(A) = -[</.(A,Ao)-0(A,±(?)] (3.3) 



*It is clear that the bare counterpart ito — Po(A) — ^£o(A) of this function fulfills this property: Mq admits a 
unique zero Ao on R, and moreover «o (Aq) < 0. Considering the form of the integral equations (|2.15p and (|2.16p . 
one easily sees that u' admits at least one zero. The case of several saddle-points can in principle be treated 
similarly, and will give rise to several contributions. 

^We do not consider here the case Aq — ±g. Indeed, this specific case would demand some further analysis. 
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between the ground state and an excited state with one particle at Aq and one hole at ±q in 
()3.ip (space-like regime), and the shift functions 



(3.4) 



between the ground state and an excited state with one particle at zbg and one hole at Aq in 
()3.2p (time-like regime). 

The corresponding amplitudes are given in terms of some properly normalized form factors 
Tjil of the current operator between the ground state and an excited state containing one 
particle at iip and one hole at fif^. The connexion 



LF%'= lim —) 



(V'L|j-(o,o)|Vc; 



(3.5) 



between the amplitude \J'lq\ and the square of the norm of an Umklapp form factor of the 
density operator between the ground state | ^lJg ) and an excited state | ^l''!_q ) containing one 
particle and one hole at the opposite ends of the Fermi boundaries was already noticed in 
[451 146| . A similar phenomenon happens for the other terms. More precisely, the amplitudes 
\j'±q\ appearing in (|3.1|) (space-like regime) correspond to the thermodynamic limit of the 
properly normalized norm squared of the form factor of the current operator taken between the 
ground state | ^pg ) and an eigenstate | tp^^ ) with a particle at Aq and a hole at ±q: 



■Ao|2 



Ao|2 



lim — 
L->+oo \27r/ 

lim — 

L^+oo \ 2ttJ 



I i(0, 0)1^^9 



ii(o,o)iv9: 



U-'gW-Ual 



(3.6) 
(3.7) 



Similarly, the amplitudes j-^^'^P appearing in (|3.2|) (time-like regime) correspond to the ther- 



modynamic limit of the properly normalized norm squared of the form factor of the current 
operator between the ground state | ipg ) and an eigenstate | ip^^ ) with a particle at ±q and a 
hole at Aq: 



I ti I 

KaoI 



KAo 



L \l+lFl(q)+l? + lFl{-g)]^ 



lim 



' L^+oo\2irJ 



{^p-;;\j{0,0)\^Pg 



The explicit expressions of p.5p -( f3^ are given in Appendix lA. 31 



(3.8) 
(3.9) 



We now comment this result. First of all, the first two lines of each expression appear as 
a direct generalization, involving the relativistic combinations x it vpt, of the large-distance 
asymptotic behavior of the corresponding static correlation function. They are in agreement 
(up to the amplitude of the third term, computed in |45| Il6] ) with the predictions coming 
from the CFT/Luttinger liquid approximations. This is indeed not surprising, since such terms 
involve particle/hole excitations on the Fermi boundary and are therefore correctly taken into 
account by a linearization of the spectrum around this point. However, when time is of the same 
order of magnitude as distance {i.e. x/t ~ 0(1)), other types of excitations also contribute to 



9 



the asymptotics, namely those involving particles (for the space-like regime) or holes (for the 
time- like regime) with rapidities located at the saddle-point Aq- These contributions appear in 
the last two lines of dSlJ, 

Remark 3.1. This result is also valid when x/t — )• (t >> x). In the t/x — )• limit {i.e. t « x), 
one has Aq — +oo, and the last two terms of (I3.ip - (l3.2p exhibit a very quick oscillation except 
around x = 0, which reconstructs the J(x)-part of the equal-time correlation function. Hence, 
as expected, the contribution of these terms vanishes at large distances. 

Remark 3.2. It is easy to see that our result coincides, at the free fermion point c = +oo, 
with the leading asymptotic behavior obtained through a direct computation of the correlation 
function. In that case Aq = x/{2t) and the combination of the last two terms of (j3.ip - (j3.2p 
reduces to: 

iJ-e'°fe*"(tF+*'^) ^ , (3.10) 

V t \ x-2qt x + 2qt y ^ ^ 

where a = 1 in the space-like regime and a = — 1 in the time- like regime. Note that this 
contribution is dominant with respect to the other power-law decaying terms. 



4 Form factor expansion for the correlation function in finite 
volume 

The derivation of the result we have just announced is based on the form factor expansion for 
the correlation function. In this section, we write down the form factor series representation 
for the current-current correlation function in finite volume, or more precisely for its generating 
function. Each term of the series can be expressed by using the finite-size determinant repre- 
sentation for the form factors of the model and their relation with the overlap scalar products. 
Such a series will be the starting point for our study. 

Let us consider the zero-temperature correlation function of currents ( j(x, t) j(0, 0) ), with 
j{x,t) = e**^j(x, 0) e~**^. Inserting the sum over the complete (see [23]) set of Bethe states 
I ^' ) between the two operators we obtain 

in which i5ex denotes the difference of energies between the excited state | ip' ) and the ground 
state \ipg). The matrix elements of j{x, 0) can be easily obtained in terms of the matrix elements 
of the operator e'^'^^, with Qx = Jq j{y,0)dy : 



{i,g\j{x,0)\i;') „ ^ (V'({A,})|e''«^|V'K(K}): 



(4.2) 



13=0 



In their turn, the latter can be computed from the results of |44| [771 158] in terms of the scalar 
products: 

(V'({A,})|e/'«MV'.(K})) =e^^^-(V'({A,}) | Mif^i,})), (4-3) 

with = '}2f=i\Po{f^ij) ~ Poi^j)]- 111 these expressions, {A} parametrizes the ground state 
whereas {fig } is a set of solutions to the twisted Bethe equations (|2.7p associated with the set 
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of integers such that | ipnifJ'ej) ) — ^ IV'') when At — )• 1 (we recall that k = e^). Note that 

it is enough to consider states with N,^ = N, as otherwise, the corresponding matrix element is 
zero. 

Therefore, (14. ip can be rewritten as 



{j(i,t)j (0,0); 



£i<...<£n 



llV'({A,})||-||V'.(K})||j,,^=o 



, (4.4) 



in which we have set Vex = limre-i-i'P^- Using the fact that Vex = if | ) = \ ) 

and that \ ipnHf^e })) = if \ 7^ \ 'ip{{fii,})) in virtue of the orthogonality 

of Bethe states [23l, we obtain 



(i(rc,i)i(0,0)) =-[85P«|„.„] 



^ ex 1 1„;. 1 12 I l„;./l 12 ' V*''-'J 



IIV'.IP-IIV' 



'112 



On the other hand, setting £1^-^. = ^^i[eo(^£a) ~ ^o(Aa)] and using similar arguments. 



2;,2 J ^ixVS^-ites^ 



did- 



IIV'«(K.})I|-||V'({A,})|| 



/3=0 



-2(l-(5^ 



1/1', i/i. 



v 



5/3(V'({Aj}) I V'K({^i;,}))/3=0 • 9/3(V'k({/^£,}) I V'({Ai}))/3=0 



IIV'.IP-IIV'^ 



,/||2 



(4.f 



Above we have used the fact that, for /3 G iM, the solutions of the k- twisted Bethe equations 
(12.71) are real, which means that 



(v^^A}) I V'^iK})) = (V'.a^^,}) I 

Let us define 



(4.7) 



£i<-<£iv 



«^E^=l["o(Wa)-"0(Aa)] 



(V'({A,})|V'.(K})) 



|V'«(K.})||-IIV'({A,})|| 



(4. 



with tio(A) = po(A) — ^eo(A). Comparing (I4.5p and (I4.6p . we get that Q^(x,t) is the generating 
function of the time and space-dependent zero-temperature current-current correlation function 
of the NLS model in finite volume: 



1 



(j(x,t)j(0,0)) =|a29|Q5v(^,i)|;3=o- 



(4.9) 



We recall that, in (14. Sp . the set of parameters {Aj} denotes the solution of the Bethe equations 
(|2.2p parametrizing the ground state of the Hamiltonian (|2.ip . whereas {/i^^} stands for the 
unique solution of the K-twisted logarithmic Bethe equations (12. 7p defined by the choice of 
integers i\ < ■ ■ ■ < In- The sum runs here through all the possible choices of ordered A^-tuples 
of integers £1 < ••• < ^^r, and therefore through all the excited states |V'K({^^j})) the 
A'^-particle sector [23j. 
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Each term of the form factor series (|4.8p involves a normahzed scalar product between a 
Bethe state and a twisted Bethe state. For the model in finite volume, such scalar products (and 
form factors) admit finite-size determinant representations |81l I58| . When the size L becomes 
large, these scalar products (and form factors) exhibit a non-trivial behavior with respect to 
the size of the system [82\ [2]. For scalar products between the ground state and an excited 
state with a finite number of particles and holes such as those described in Section [21 it was 
shown in |46l H3] how to extract the leading large L asymptotic behavior from their finite size 
determinant representations. 

The detailed representation for the products of form factor that appear in each term of the 
series M.Sh . as well as its leading behavior in the thermodynamic limit, is recalled in Appendix lAl 
In finite volume, the corresponding normalized scalar product takes the following form: 



(v({A,}) I v.(K}): 



|V'.(K})||-||V'({A,})|| 



^Jv({Ai},{/^£,}) • ^Jv({Aj},{At£,}), 



in which 



N 



DN{{Xj},{f^i,}) = ll 



sin2[7rF(Aj 



det 



N 



A,- 



^4 



(4.10) 



(4.11) 



is the so-called discrete part of the form factor, which contains the whole non-trivial singu- 
lar part of the form factor with respect to the system-size (and is quite universal), whereas 
^Af({Aj}, {f^ej}) is a dressing function which, for a n-particle/n-hole excited state as defined in 
Section [21 admits a smooth thermodynamic limit: 



lim gN{{Xj},{fJ'e,}) 



(4.12) 



Here Qn is a holomorphic function of the particle/hole rapidities and {fJ^ha}- 

croscopic details of this dressing function (see Appendix [A]) depend on the model and do not 
really matter for our study. We will only use the fact that Qn can be represented in terms of a 
functional 



{^J'Pa} 



of the function 



X 



{^J'ha} 



E 

a=l 



A 



A 



(4.13) 



(4.14) 



The behavior of the thermodynamic limit of the discrete part is more difficult to obtain 
(see [461 '43j and Appendix [X]). Its knowledge is however unnecessary for our purpose, since our 
study will directly rely on the finite-size formula (14. lip . We will simply use the fact that (14. lip 
can be understood as a functional of the finite-size shift function F, depending moreover on the 
extra sets of particle/hole integers {pa} and {ha}'- 



-CAf({Aj},{/ifJ) = Dn,, 



{Pa} 
{ha} 



(4.15) 



More precisely, the function ^ being fixed and the set of parameters Xj being defined as the 
pre-image of the set {j/L : j = 1, . . . , N} by this function, it means that: 
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• the function in (|4.11|) should be understood as a functional of the shift function F 
through the relation = — L~^F] 

• the set of integers ij being fixed by the choice of {pa} and {ha} (and vice-versa), the 
parameters fi^. are obtained as their pre-image by the function L^^- 

5 Effective decoupling of the finite-size form factor series 

The structure of the summations in (|4.8p is extremely intricate. Indeed, the sets {ni-} of 
Bethe parameters over which we sum up are implicit functions of the N integers ij labelling 
the corresponding excited states. Although one can build on such a description so as to use 
the multidimensional residue theory to recast this complex sum into a contour integral (the so- 
called master equation later still remains difficult to handle. In particular, one 
cannot perform the thermodynamic limit directly on the level of the master equation. Therefore, 
one has to expand the contour integral into yet another series (the so-called multidimensional 
Fredholm series) that has a presumably well defined thermodynamic limit. The asymptotic 
analysis of the obtained series, decomposed in terms of cycle integrals, then relies on a non- 
trivial summation process, resulting in particular into the dressing of bare quantities (energy, 
momentum) into dressed ones. Although such an approach was successfully applied in [45j to 
produce the leading asymptotic behavior of the time-independent correlation functions of the 
XXZ chain and of the NLS model (see also [63] for the temperature-dependent case), we wish 
here, so as to bypass some of the subtleties of |45| . to propose an alternative line of though by 
pursuing our study directly on the form factor series. 

The idea of our method is the following: the series (j4.8p can be related to an essentially 
free fermion type series (see Appendix [B]), this via some reasonable physical assumptions con- 
cerning the contribution of each of its terms in the thermodynamic limit, and via some formal 
manipulations. Indeed, it is shown in Appendix |B] that if, for each excited state, 

• the rapidities of the particles and holes are decoupled, i.e. they are determined by a 
counting function that does not depend on the position of the roots of the corresponding 
excited state, 

• the dressing part ^at of the form factors is decoupled, 

then the series in finite volume can be summed up into a finite-size determinant. We explain in 
this section how to reduce the study of our highly coupled series to this simple case. It means 
that, in our setting, the key-role in the summation is played by the universal discrete part Dj^ 
of the form factors exactly as in the free fermion case. This fact should be put in parallel 
with the key- role played, in the master equation-based asymptotic analysis |45| . by the Cauchy 
determinant part of the form factors. 

5.1 An effective form factor series 

Considering that we will finally be interested in the thermodynamic limit L — )• -|-oo of the 
generating function (j4.8p . we now slightly simplify its form factor series. Our simplifications 
rely on two kinds of assumptions: 

(i) the contribution to the sum (|4.8p of a state having a macroscopically (with respect to L) 
different energy and momentum from the ground state should vanish in the thermody- 
namic limit; 
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(ii) for each term of the series, "non-discrete" quantities (i.e. parts of the form factor expres- 
sion behaving smoothly at the thermodynamic hmit) should contribute to the thermody- 
namic limit of the sum (j4.8p only through their leading order in L. 

Assumption (i) can be attributed in particular to the extremely quick oscillation of the phase 
factors for states having large excitation momenta and energies. It means notably that: 

• we can use the particle/hole picture to describe the large-L behavior of the form factors, 
keeping in mind that the relevant (i.e. the one not vanishing in the L — )• -|-oo limit) part 
of the form factor expansion corresponds to a summation over states with a finite (or at 
least growing much less than L) number of particle/hole excitations above the ground 
state; 

• we can introduce a "cut-off" (with respect to L) of the range of integers on which we sum 
up in (j4.8p since, for very large integers, momentum and energy of the corresponding state 
become also very large. 

Doing this, we roughly speaking neglect correcting terms in the lattice size L. It is thus reason- 
able to assume that, on the same ground, one can also neglect some finite-size corrections to 
the leading thermodynamic behavior of the form factors, as stated in (ii). More precisely, we 
make the following simplifications in the series (14. 8p : 

• we replace the finite-size shift function F, the smooth part Gn and the oscillating exponent 
^a=li'^o(^^ia) ~ '"o('^a)] by their leading contribution in the thermodynamic limit; 

• in the obtained expression, we now understand the rapidities of the particles fip^ (resp. 
of the holes ///i^), corresponding to a particular choice of integers ii < ■ ■ ■ < In, as being 
defined in terms of the thermodynamic limit F of the shift function (and no longer in 
terms of its finite-size counterpart (12. 9p ) as the pre-images oipa/L (resp. ha/L) by the 
counting function = £, — F/L; in other words, the 2n parameters {/Up^} and {/x^^} are 
obtained as the unique solutions to the system of 2n equations 



Pa , 

— and 
Ij 



with a = 1, 



, n. 



These simplifications result into an effective series that we assume to have the same value 
at the thermodynamic limit as the original series. In other words, we conjecture that. 



lim Q^(x,t)efr 

L,N-^oo 



lim Q%{x,t) = Q''{x,t), 



(5.2) 



with 



E 



h<-<tN 



F 



X ^ 



il^ha} 



Wr, 



{l^Pa) 
{ll'ha} 



. (5.3) 



In (15. 3p . the sums are restricted to the set ;Bl = G Z | —wl < ^ < ^^l}, where wl is 
a cut-off growing as L^^'', for some e > 0. For a given set of integers {^j}, the parameters 
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Pa and ha define the position of particles and holes as explained in Section [21 The (finite- 
size) representation for the form factor has been partly replaced by its thermodynamic limit. 
More precisely, the dressing function and the phase factor, which both behave smoothly at 
the thermodynamic limit, have been replaced by their leading equivalents. Note that, as a 
consequence, it is now the dressed energy and momentum that appear in the phase factor (c/ 
Appendix [a]): u{X) = p{X) — ^^(A). In (15. 3p . we have kept the finite-size expression of the 
discrete part, but we have nevertheless replaced the expression of the finite-size shift function F 
by its limiting value F ()2.14p . Finally, the particles and holes' rapidities associated to a given 
set of integers ii < ■ ■ ■ < £n are obtained (in terms of F) via the system of 2n equations (15. ip . 

Remark 5.1. We stress that, in (15. 3p . we are still in "finite volume" {i.e. for the moment, L 
and are kept finite). We have only modified the original series (j4.8p arguing that it should 
not affect the value of its thermodynamic limit. 

Remark 5.2. The introduction of the cut-off wl is convenient to avoid problems of convergence 
in our further manipulations of the series. Indeed, as far as we remain in finite volume, we now 
deal with finite sums only. 



5.2 Towards a free-fermion type series 



Although the effective series (|5.3p is already simpler that the original one it is still highly 
coupled: 



• the thermodynamic limit of the shift function (|2.14p still depends on the position of 
particles and holes, 

• the rapidities of the particles and holes have to be computed for each excited state sepa- 
rately (i.e. they are excited state dependent), 

• the expression of the functional 1^ is extremely intricate. 

Our aim is now to relate QAr(^)^)cfr ^° ^ decoupled series, so as to reduce its analysis to the 
one of the generalized free-fermion case that is carried out in Appendix [Bj This can be done 
by understanding the coupling between variables as the result of the action of some functional 
translation operator^ (see Appendix [U]) . 

The function tz7„ (|4.14p depends linearly on a function of the rapidities of the particles and 
of the holes. This means that one can formally express this dependence by means of the action 
of a functional translation operator, as explained in Appendix O for any sufficiently smooth 
functional J- supported on a neighborhood U of the real axis. 



J" 



{iJ'ha] 



dA 



a=l 



1 



1 



5uj{\) 



(5.4) 



1.^=0 



The thermodynamic limit of the shift function (I2.14p also depends linearly on a function of the 
rapidities of the particles and holes. However, the situation is slightly more complicated than in 
(j5.4p since, in virtue of (jS.ip . the parameters {/ip„} and {^J-ha) ^'^^ themselves functionals of the 
shift function. As explained in Appendix [Cl one can still represent any smooth functional of 



®The use of functional translation operators is simply a convenient and compact way to manipulate generalized 
Lagrange series (see Appendix C of [45]). In fact, the whole reasoning that follows can instead be performed by 
writing explicitly the corresponding series. 
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this shift function in terms of translation operators, provided that one imposes some operator 
ordering : • : and that one takes into account the contribution of the Jacobian coming from 
the summation of the corresponding multi-dimensional Lagrange series (see formulae (|C.5p and 
(IC3ID): 



T 



F 




dA [0(A,^ftJ - (/>(A,AipJ] 



6t{\) 



J. (5.5) 



Above the operator ordering : • : is such that, in the formal series expansion in powers of 
(5/(5r(A), all functional derivative operators are located on the left, v^- is the function 



-.(A)=^^+r(A), 



(5.6) 



and J is the Jacobian J = detiR [/ — 5T[ur\{\) / 5t{^)\^ w-f °^ functional 



rK](A) = 5^[</.(A,/z,j 

a=l 



'/'(A>A*Pa)]> where iii> = l, 



-1 



L 



with 



(5.7) 



This Jacobian is evaluated at r = -F, where F coincides with the shift function occuring in the 
l.h.s. of ()5.5p . We stress that the bar occuring in the parameters Jip^ (resp. Jlha) appearing in 
the r.h.s. of (j5.5p indicates that these are to be understood as the pre-images of Pa/L (resp. 
ha/L) by the counting function ^y^, as in (15 .jp . 

It is easy to convince oneself that, provided that the number of particle/hole excitations is 
finite, J = 1 + 0(-L~^). In the light of our previous arguments, only this sector of excitations is 
expected to contribute to the L — )• +oo limit of the form factor expansion. As a consequence, in 
order to be consistent with our previous approximations, we also drop the finite-size corrections 
due to this Jacobian by considering from now on that our effective series (still abusively denoted 
by the same symbol Q^(2;,t)eff) is in fact given by the expression (15. 3p in which each term is 
multiplied by the inverse of this Jacobian. 

Hence, the use of the functional derivative leads to the following representation: 



DN,n 



{Pa} 
{ha} 



F 



J 



-1 



;llgE Me* 



a=l 



a=l 



M-)l : 



T=0 

aj=0 



. (5.8) 



Here we have insisted on the fact that D]\j n and $f are functionals acting on the argument • 
of ^'r(-) and ti;(-). Also, on the l.h.s. of (j5.8p . we have explicitly pointed out the parametric 
dependence of F and tZ7„ on the parameters {/ip„} and 

In the r.h.s. of (15. Sp . the functional D^^^n acts on a shift function that does not depend 
anymore on the particle/hole rapidities. In other words, the effective shift function v-r becomes 
independent of the summation variables and hence mimics the one appearing in the generalized 
free fermion model studied in Appendix |Bj Moreover, the equations defining the position of the 
particle's/hole's rapidities also become decoupled: the rapidity Ji^ is now the unique solution 



16 



to (^u-j-ifii) = ^/L, and its value does not depend anymore on the choice of the other integers 
describing the excited state, but only on the function r. 

It is convenient to express, on a formal level, each term of the pre-factor in (jS.Sp as a ratio 
of two exponents ediP-pa) / Qdifiha) ^ where g is an operator valued function (we have used the hat 
so as to insist on this property): 



9 = 9T + guj, with gr{X) 



dfi 5 
fi — X 6uj{X) 



(5.9) 



In the following, we set 



-ixu—g 



(5.10) 



The operator order : • : being linear, we can exchance it with a finite summation symbol 
such as the one appearing in the expression (j5.3p for Qpf{x,t)gfi. Hence, we obtain 



I3xpp 



a=l E-iP-Pa 



D 



N,n 



ei<-<iN 



{Pa} 
{ha} 



T=0 

a)=0 



a=l ^-\Xa) 



T=0 

aj=0 



(5.11) 



in which X^r is the generalized free-fermionic functional (jB.ip studied in Appendix [Bl 



5.3 Representation in terms of a finite-size determinant 

In Appendix [Bl we have recast the generalized free-fermionic functional Xjv[j^, -E^] into a finite- 
size determinant ()B.2p . this without any further approximations. Such a representation is 
obtained through "purely algebraic manipulations on the initial definition (IB.ip for Xjv[i^, -E^]. 
These two representations are still equal when applied to the (non-commutative) operator valued 
functions Vt and -E_ : it is indeed not a problem to implement the appropriate operator order at 
any step of the computation performed in Appendix [B] One may therefore use the determinant 
representation (IB.2P to obtain new representations for our series (15 . 1 1 p . provided that one 
imposes an operator order on the entries of the columns of the determinanl|3 [84]. This leads to 



Q^(a;,t)cfr 



N 



5jk + 



H'{Xk) 



T=0 

uj=0 



where the (operator valued) kernel is given as 



V^''\Xk,X,) 



sm[iTu-r{Xk)] sin[7ri/^(Aj)] r P,(L) 



2i7r(Afc — Xj 



E\^>iXk)E^{X, 



Ei^\x,)E^{Xk)}, 



(5.12) 



(5.13) 



with 



E^_^\X) = iE_{X) 



— — — y + COt[7rUr{X)\+ lL[l^r,E_ \{X) } . 



(5.14) 



^Such an order can for instance be implemented explicitly by expanding the determinant into a sum over 
permutations. 
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The expression of the functional II can be found in (jB.Sp . 

Such an expression provides a formal re-summation of the effective form factor series (I5.3p . 
In the next section, we show how to take the thermodynamic limit of (j5.12p . 



6 Thermodynamic limit and asymptotic analysis 

We now build on the results of the previous section so as to derive the leading large-distance/long- 
time asymptotic behavior of the thermodynamic limit Q'^{x,t) of the generating function 

Q^(x,t). 

6.1 Representation in terms of a Predholm determinant in the thermody- 
namic hmit 

The next step of our study is to obtain a convenient representation for the thermodynamic limit 
N,L ^ +00 of the generating function (|4.8p . Recall at this point that the effective series (|5.3p 
was built so as to approach, in this limit, the same value Q'^{x,t) as the original form factor 
series (|4.8p . We will therefore use the representation (j5.12p to derive a suitable expression for 

Q"(x,t). 

In principle, prior to taking the thermodynamic limit of a formal representation such as 
(I5.12p . one should first compute the effect of the translation operators: the operation of taking 
the thermodynamic limit is indeed a priori only allowed on expressions defined in terms of 
explicit holomorphic functions (i.e. which do not contain any operator valued functions). In 
Appendix ini such a procedure is explicitly performed so as to obtain, starting from (|5.12p . a 
particular series representation of the thermodynamic limit Q'^{x, t) of Q'^{x, t). More precisely, 
in Appendix ini we expand the finite-size determinant in (|5.12p into a sum over determinants of 
all sub-matrices of V{Xj, A^), factorize the Cauchy part of each of these determinants, compute 
the effects of the translation operators, and then take the thermodynamic limit. This enables 
us to express Q'^{x,t) as a series of multiple integrals coinciding, in the equal-time case, with 
the series expansion derived from the master equation representation in |45) . 

In fact, it is easy to convince oneself, by carrying out the computation in the reverse order, 
that the operation of taking the thermodynamic limit A^, L — )• +00 and the one of computing the 
effect of the translation operators do actually commute. Indeed, let us consider the expression 
(I5.12P in which we have sent directly N,L ^ +00, hence replacing sums by integrals in the 
prefactor and the finite-size determinant representation for X]\f by the Fredholm determinant 
([RTT]1 : 



-p.pp - J Ur{X)[ixu'{X)+g'{X)]dX 

: e det [I + V[iyr,u,g\\ -Wluj] 



5.1) 



where 



with 



2iir{X — jj) 

■ 9(A) 

E_(A) =e-*^ 2 - 2 , (6.3) 
E^iX) = ^E_{X) { / + ^4^cot[..(A)] ) . (6.4) 
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Then, if we decompose the Fredholm determinant in (|6.ip into its Fredhohn series and subse- 
quently compute the effect of the functional translation operators, we obtain the same represen- 
tation (|D.24p as by performing all these operations in the opposite order. This means (provided 
that the series (ID.24P is convergent, which is not completely obvious but seems nevertheless 
a reasonable assumption) that the thermodynamic limit of the effective series (j5.12p . which is 
supposed to give the thermodynamic limit of the original form factor series (14. 8p , is effectively 
given by the expression (|6.ip . i.e. 

_0 - f Ur{X)[ixu'iX)+g'{\)]d\ 

Q^{x,t)=e^ : e det[l + V[u^,u,g]]-W[u;]:^^^, (6.5) 

uj=0 

with a kernel V given by (|6.2p -( l6l^ . It also means that we can now use any other existing 
representation for the Fredholm determinant so as to compute the effect of the translation 
operators and recover standard scalar- valued functions. In fact, it is not very convenient for our 
purpose to expand the determinant into its Fredholm series like in Appendix |D] since the latter 
does not provide any information on its large- asymptotic behavior. This large- j; asymptotic 
behavior was studied in ||60j precisely with the goal of computing the asymptotic behavior of 
()6.5p . and a much more convenient (with respect to the x — )■ -|-oo limit) representation for the 
Fredholm determinant was obtained there. 

6.2 Large-x asymptotic behavior of the Fredholm determinant 



The large x asymptotic analysis of the Fredholm determinant with kernel (I6.2p - (l6.4p was per- 
formed in |60| using Riemann-Hilbert problem-based techniques. There, it was proven that, 
under some hypothesis about the regularity and behavior of the functions u, u and g defined on 
some open neighborhood of the real axis (see [60j for more precisions), and provided that the 
function u has a unique saddle-point Aq on R (with Aq / ±g), 

det[_g,g] [/ + V] [u, u, g] = exp J y [ixn'(A) + ^'(A)] i/(A) dA I 



[ e=±l 



v + e,u\ 



+ \Y^ e^"^['^(^«)-"(^«)]+°[s(^«)-^'(^5)l6S''"^[i/,u]e^[zy,u] +7^:,.[^^,M,^ (6.6) 

e=±l J 

with a = +1 in the space-like regim^ (Aq > q) and a = —1 in the time-like regime (Aq S] — q\). 
In (16. 6p . the functional Bx reads 



*The case Ao < — g was not considered in [60) . 
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M)r., ..1 [2gxn'(6<7)+^0+]2--M {-e)i^{eq) ( Aq + (? \-2-(Ao) 



with given by ([02]) . and 

[j, yl = e"' 

^ V-2vrn"(Ao)u'(eg) {\o - eqf \\a - q - m+ , 

( -2inuiXo) _ 1 xl-a / -2i7ri.(£g) _ -, ^a r(l - eaiyjeq)) aJlu]iXo)-aJlu](eg) on 

^ ^ ^ ^ r(l + eaKe'7)) ' ^ ' 

where J is given by ()A.25p . TZxW^u^g] is a remainder which is uniformly of order 0(^^^) in 
what concerns the non-oscillating corrections, of order 0(i3a;[z^ + e,u\ ^^^) in what concerns 
the oscillating corrections at q^'^^1-^(i)~'^(~i)^ ^ and of order 0(^Bx[i^,u\b^i''^\i',u\^^§jj^ in what 
concerns the oscillating corrections at eiaxHXo)-~u(eq)] ^ 

was shown in [60j how to obtain a 
series representation for this remainder, the so-called Natte series, which possesses the property 
of being well ordered with respect to its large-x behavior (in can be shown that its n-th term 
is (uniformly) a 0(x~"'*), for some < a < 1 that depends on v, u and g). Hence, this series is 
well adapted for the study of the asymptotic behavior of the determinant (on the contrary to 
its Fredholm series). The form of such a Natte series is recalled in Appendix [El 

Remark 6.1. Note that for |K {u)\ < 1/2 the first (non-oscillating) term Bxy^, u] in (16. 6p is always 
leading, at large x, with respect to the other ones. This leading term will produce the leading 
asymptotic behavior of the generating function Q'*(x,t), as we will see in the next subsection. 
However, recall that we have to differentiate twice with respect to x and with respect to (3 at 
/3 = in order to obtain the correlation function ( j(x, i) j(0, 0) ). By such a process, the first 
oscillating corrections may become leading with respect to non-oscillating terms. This is the 
reason why we also consider these corrections in (|6.6p . 



6.3 Asymptotic behavior of the correlation function 

So as to obtain the asymptotic behavior of the correlation function, it remains to compute the 
effect of the functional translation operators on the representation (|6.6p of (|6.5p . First, we 
observe that the exponential pre- factor of ()6.5p is canceled by the one in (I6.6p : 



exp I - j Urifi) g'in) dfj. > exp | j zy^(/i) ^(^) d^ i = 1. (6.9) 



-q --q 



Such a simplification is justified by the fact that, when applying (16. 9p on any functionals of 1/^ 
or u, we observe an algebraic cancellation for each term of the : • : ordered series expansion in 
5/5t and S/duj for these exponents Therefore, 

I e=±l 



e=±l 



r=0 
u)=0 



(6.10) 



*Such an algebraic cancellation is very similar to the one occuring when computing — A)'^ through a 
Taylor series expansion around zero. 
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It can be proved that, by computing the effect of tlie functional translation operators oc- 
curing in the remainder TZxli^TjUjg], one obtains corrections that are of the same order {i.e. 
0(^2££) corrections to each of the terms already present in the asymptotics) as originally in 
()6.6p . This is explicitly done in Appendix IeI by using the so-called Natte series representation 
of the Fredholm determinant derived in [60j. More precisely, it is shown in Appendix [El that 
the effect of the translation operators do not mix the orders in x among the different terms of 
this (well-ordered) series. Therefore, the leading asymptotic behavior of the generating function 
follows directly from the above leading asymptotic behavior of the Fredholm determinant. 

Since no translation operator is applied on the first (non-oscillating) term of (j6.10p . we 
simply need to set r = and a; = into the corresponding expressions. The action of the 
translation operators e''l9i<i)-9(-<i)] on the second term of (I6.10p results into replacing the function 
^[w] by 



-eq 



-eq 



and the function Ur by the shift function 
eq \ ifi 



F 



-eq 



27r 



-Z 



[<^(-,6(?)-0(-,-6.?)] = (^+6)z 



(6.11) 



(6.12) 



associated to a state with one particle and one hole located on the opposite ends of the Fermi 
zone. Similarly, the action of the translation operators Q°'\9{^o)-g{<^q)] results into replacing '^[uS\ 
and v-r respectively by 





AoV 







and F 



in the space-like regime a = +1, and by 



= Q 

in the time-like regime a 



and F 



iJ3_ 
27r 



i/3 



Z -[(/>(•, Ao) -(/.(•, eg)] 



^Z-[c/>(.,6g)-0(-,Ao)] 



(6.13) 



(6.14) 



-1. Therefore, we get 



Q-{x,t)=e-^{ Bx 



— Z, u 
2-K ' 



e=±l 



^0 ( 1 + O 



he Z, n 

27r ' ' 



F 



log X 

X 

eg 
-eq 

F 



1 + 



logx 



1 + 



log X 



(6.15) 



in the space-like regime, and similar expressions in the time-like regime. Namely, in the time- 



like regime, the exponent in the last term changes sign and h^^'^^^ , " | eg ) ) (eg) ^"^^ replaced, 

respectively, by h^{^-^\ F( • |fj and 

In order to obtain the leading asymptotic behavior of the correlation function of currents, 
it remains finally to compute the second x-derivative and second /3-derivative at /3 = of the 
previous result, cf (j4.9p . The derivatives of the first term (the non-oscillating one) in (|6.15p 
produce the constant and the non-oscillating term appearing in (|3.ip and (|3.2p (we have used 
(IA.20P ). In their turns, the derivatives of the two types of oscillating corrections produce the 
corresponding oscillating terms in (|3.ip . (|3.2p . with amplitudes given by (|A.28p - (|A.30p . 
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7 Conclusion 



In this article, we have proposed a new method to derive, starting from first principles, the 
leading asymptotic behavior of the two-point correlation functions of quantum integrable sys- 
tems. To explain the main steps of this method, we chose to focus on the example of the 
current-current correlation function of the quantum non-linear Schrodinger model. The case 
of the field/conjugated field correlator ( t) ^'(0, 0) ), together with a rigorous setting for 
carrying out all the manipulations with operator valued determinants, will be given in [59J. 

Our result goes beyond the CFT/Luttinger liquid based predictions: the saddle-point con- 
tributions that appear in the asymptotic behavior (|3.ip - (j3.2p involve excitations away from the 
Fermi surface and cannot be neglected for xjt finite. 

Compared to the approach of |45| . based on the master equation representation for the 
correlation functions, the present study relies directly on their form factor expansion. We 
would like to conclude this article by making here a few comments about the similarities and 
differences between the two methods. 

Of course, since the master equation can be understood as the result of a summation over 
the form factors, the spirit of the two approaches is essentially the same. In particular, in 
both approaches, the asymptotic analysis relies essentially on the singular part of the form 
factor which is explicitly extracted (in the form of a Cauchy determinant squared) , whereas the 
(mo del- dependent) regular part is treated as a dressing part which is formally decoupled (or 
linked to decoupled functions) to make the analysis possible. This enables one to draw a link 
between the quantity to estimate and the Fredholm determinant of a generalized sine kernel. 
The asymptotic analysis of this Fredholm determinant then leads to the asymptotic behavior 
of the correlation function. 

However, there exist some essential differences between the two approaches. In |45| . the 
correlation function was expanded into a series whose building blocks were the so-called cycle 
integrals. These cycle integrals could then be related to a Fredholm determinant, which allowed 
one to access to their asymptotic behavior. The physical interpretation of these objects was 
however not clear, and they happened to be quite indirectly related to the correlation functions. 
In fact, once the asymptotic behavior of these cycle integrals established, one had to sum up 
the series so as to obtain the asymptotic behavior of the correlation function. In this process, 
the main problem was that not only the leading asymptotic behavior of individual cycle integrals 
was contributing to the leading order for the correlation functions: one had to perform a fine 
and non-trivial study |47| SH [62] of the asymptotic series so as to gather terms at all order 
that finally, when summed up, were contributing to the leading order. In fact, all these terms 
were rearranging themselves into some generalization of a multiple Lagrange series producing, 
when summed up, a dressing of bare quantities (energy and momentum) into dressed ones. 
On the contrary, here, one deals from the very beginning with objects having a clear physical 
interpretation, the form factors. In this context, the dressed quantities appear naturally when 
considering the thermodynamic limit of these form factors. Hence, performing the summation 
over the form factors, we can connect the series to a Fredholm determinant that is already 
expressed in terms of these dressed quantities. Therefore, the asymptotic study is much simpler: 
it happens that the leading asymptotic behavior of the Fredholm determinant gives directly the 
leading asymptotic behavior of the correlation function. In other words, there is no need to 
resort to highly non-trivial summation as in |45) . 

In fact, we would like to stress that the whole process described in the core of this article 
is quite simple and direct. Once the form factor series written down and the effective con- 
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tributing part of the form factor established, the result of the summation can be expressed, 
in the thermodynamic limit, in terms of a Fredholm determinant. The leading asymptotic be- 
havior of the correlation function follows directly from the leading asymptotic behavior of this 
Fredhom determinant. Even if, in the course of the computations, we use some functional trans- 
lation operators to relate our series to a decoupled one, whenever the action of these functional 
translation operators has to be computed, it is quite straightforward, i.e. it produces a simple 
translation of the functional on which it acts. On the contrary, if one wants to recover, starting 
from the form factor expansion, a series similar to the one studied in |45| . the computations are 
much more involved (see Appendix [P]) . In particular, the action of the translation operators 
on the series of Appendix [D] produces non-trivial effects: one has to deal with summations 
of generalized Lagrange series, which results into an undressing of the dressed quantities into 
bare ones. It is therefore clear that most of the mathematical complexity corresponding to this 
non-trivial summation has already been taken into account by the fact that we had some pre- 
cise description (the particle-hole picture) of our form factors, which allowed us from the very 
beginning to deal with dressed quantities instead of bare ones. Note that the simpler setting of 
our method enables us here to consider the time-dependent case, which is not so obvious within 
the approach of [45]. Note also that, in principle, there is no intrinsic obstruction preventing 
us from obtaining higher order terms in the asymptotic expansion for the correlation functions. 
For this it is enough to refine the asymptotic expansion of the Fredholm determinant. 

Of course, there is a price to pay for this simplicity: the need of some clear picture to describe 
the form factors. Whereas in jlS] all kinds of Bethe roots were automatically taken into account 
within the master equation framework (and this without any precise study of the spectrum), 
here we strongly rely on the fact that the spectrum of the model we consider is particularly 
simple: all excited states can be described in terms of particles and holes. In order to apply 
our method to the case of the XXZ spin chain, for example, one has also to take into account 
the contribution of complex solutions. Although it seems that, for the time-independent case, 
these solutions do not contribute to the leading asymptotic behavior of correlation functions 
(see [42j), the question remains open in the time-dependent case. This will be the subject of a 
further study. 
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A The form factors and their thermodynamic hmit 

In this appendix, we present the explicit expression and the leading thermodynamic behavior 
of the combinations of finite- volume form factors which appear in each term of the series (14. ip , 
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i.e., with the notations of Section HI of 



m\'-\\^ 



'112 



1 



2^2 } JxVg,-it£^^ 



■ (A.l) 



/3=0 



These results being the complete analogues, in the case of the NLS model, of those derived in 
for the XXZ chain, we skip the details of the computations. 



A.l Thermodynamic limit of the space and time-dependent phase factor 

The thermodynamic limit of e^^'Pcx~^^^ox generates dressed quantities. This limit being com- 
pletely smooth, it is easy to see that 

e'"^--'*^- = exp [ix [ n'o(A) F(A) dX + ixJ2 " Mf^hJ]]{l + 0(1/L)) 



a=l 



exp|-^| n^(A)Z(A)dA + i:Ef;K;UpJ-^z(/.,J]|(l + 0(l/L)), 



in which we have used the explicit form (I2.14p of the shift function F, and where u is the 
following combination of dressed quantities u{X) = p{X) — te{X)/x. This identification was 
possible due to 



u{X)=uo{X)- / no(/i)(/)(^,A)d^. 



(A.2) 



Note that, due to the fact that s'q is an odd function whereas the dressed charge Z is even, one 
has 



J dAZ(A)uo(A) = 2pF. 



(A.3) 



-q 



A.2 Representation of the normahzed scalar product 

Let {Aj} be the solution of the system of Bethe equations (12. 2p parametrizing the ground state 
of (12. ip and let {fJ-i-} be a set of Bethe roots of (|2.7p parametrizing an excited state above the 
K- twisted ground state in the A^-particle sector. Then the normalized modulus squared of the 
corresponding overlap scalar product in finite volume can be represented as 

m{^j})\M{f^eA)) ' 



^^({Ai},{M^,}) • Wjv({A,},{/i,J) •^;v({A,},{/.,J), (A.4) 



where D]\f, Wn and An are given by 

DN{{Xj},{f^e,}) - ii 2r2?/n ^T>f 



n 2 



det 



N 



Xj - nek 



(A.5) 
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det 



N 



i^(A,-Afc) 



6ii ^ 



We recall that and ^ denote respectively the excited state counting function (12. 8p and the 
ground state counting function at k = 1, whereas F is the finite-size shift function (I2.9p . Here 



9i and ^2 are some arbitrary real parameters. The integral operator I + ^Ug^^ acts on the 
closed contour T surrounding the ground state roots {Xj} and no other singularity of the kernel. 
This kernel reads 

with 

= ^- (A.9) 

CJ + «c u — IC 

Remark A.l. Although each individual term, in (|A.7p . does depend on the set of auxiliary 
parameters 9k, the overall combination An does not. This was proven in |45| . 

The factor DnHXj}, {fJ-ij}) (lA.SP is the so-called discrete part of the form factor. It contains 
all the non-trivial "singular part" of the form factor (see [43|). On the contrary, the factor 

§N{{Xj},{l^i,}) = ^N{{Xj},{f^i^})-AN{{Xj},{f^i^}) (A.IO) 

admits a smooth thermodynamic limit. It can be thought of as a dressing function. It is equal 
to 1 at the free-fermion point. 

In the case of an excited state IV'^CI^^j})) with a finite number of particles and holes as 
described in Section [2l it is easy to compute the thermodynamic limit of the dressing function 
Gn (see |43j): 



with 



n 

-2in E E {clF]{i,p,+iec)-ClF]{i,h,+iec)}+Co[F] 



n 

-2in E E {C[F]C .„ „ 

X e "-i^-i , (A.12) 



^-2^.{C[F](e,+.c)+C[F](e,+.c)} n (g^ _ + ,c)(g^ _ + 



^ ^ (1 - 6-2-^(^1) )(1 - e-2-^(^2)) (01 - f,p^ + zc)(02 - Mp, + ic) 

det [/+2i^^eJ det [l+i^Uo,] 

det^[/-^K] • (^-^'^ 
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In these expressions, F denotes the shift function F (j2.11|) . C[z^] is the rational Cauchy transform 
of the function v on [—(?,(?], 



CM(A) 



d/i u{fi) 



2m /i — A' 
whereas Cq is the fohowing functionah 



(A.14) 



dA d/^ 



(A — ^ — ic) 



2 ■ 



(A.15) 



The integral kernel Uq takes the form 



f. {uj- fipjju - fif,^ + ic) e^--^[^l(^){i^^(^ _ uj') - K^je - oj')} 



(A.16) 



Remark A. 2. The dependence on the particle/hole rapidities of Qn is twofold: explicitly in the 
above expressions, and also in the shift function F. The latter is a holomorphic function of 
{l^Pa} ^'^d {fiha}j (|2-14p . Therefore, Qn is itself a holomorphic function of {/Xp^} and {^h^}- 
This function can, in fact, be understood as a functional ^[wn] of the function Wn (I4.14p (note 
that the dependence in n is exclusively contained in Wn)- For this it is enough to observe that 



F(A) = ^Z(A) 



d/i 



<l){\,ll)'Wn \ fl 



(A.17) 



where the contour FfM) surrounds the real axis counterclockwise, and that 

dA 



n 

a=l 



,/(Mpa)-/(A'ha) 



exp 



2m 



f{X)Wn[ A 




for any holomorphic functions / with a sufficiently mild growth at infinity on r( 
we have 



n 



dA d/i 



exp 



{2mf 



- ic) 

log(A - n-ic) Wn 



(A.18) 
In particular. 



{/^Pa} 



VJr, 




(A.19) 



Remark A. 3. The explicit factor (1 



in (IA.13P stresses that, for generic configurations of 



the parameters {/ip„} and {/U/^^}, Qn, considered as a function of /3, has in general a zero of order 
2 at /3 = 0. Note however that there exist some exceptions, in particular at the free fermion 
point c = +00, when Q is identically 1 or, for general c, when the excited state coincides with 
the ground state in the limit /3 — )• {i.e. for n = 0). In the latter case, the limit can be taken 
as in [45j and we obtain 



lim Qq 



1. 



(A.20) 
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The study of the thermodynamic behavior of the factor D^v (|A.5p is shghtly more technicah 
Proceeding as in [43j, one obtains that, for a n particle/hole excited state, 



^iv({A,},W})= Uz 



L \ 



27r J 



-FHq)-FH-q) 



B[F] 



2 f{Pj} , {Pj -N + Fi^ip.)} , {hj + FifihM ,{N + l-hj- Fifih,)} 



{Pj - N} , {p, + F(^pJ} , {hj} ,{N + l- h,} 

n „:„2r 



X det^ 



1 



sin^[vrFK)]e"[^'(^-^)-" [^'K) 



1 + 



L 



(A.21) 



The expression ()A.2ip is given in terms of the below functional of the shift function (|2.14p : 



(27r)''('?)"^("9) 
where G is the Barnes G function and 

c.M = l/dAd,'-'W''W-.^'(MHA) 

2 J A — /i 

-9 



(A.22) 



q-X 



K-g)-KA) 
q + X 



dX. (A.23) 



We have also defined 

J[F](A) = 2F(A) log 

with 



A - qp{X) -p{-q) 
X + q p{X) -p{q) 



JM(A) 
Finally, 



u{li) - u{X) 
— X 



fr r(afc) 



d/i. 



(A.24) 



(A.25) 



(A.26) 



with the prescription that, should some of the particles in (IA.2ip have their rapidities to the 
left of the Fermi zone, the arguments of the F-functions have to be understood as limits 



T{pk - N) 



lim 



r(pfc + e) 
5 T{pk -N + e) 



(A.27) 



Note that in (IA.2ip the thermodynamic limit has only been taken partly. Indeed, since the 
complete thermodynamic behavior of ()A.2ip depends on whether the corresponding particles 
and holes remain or not at a finite distance from the Fermi boundaries. In the next subsection, 
we particularize these expressions to the special form factors with one particle and one hole 
which appear in the results (|3.ip . (|3.2p . 
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A. 3 Explicit value of the amplitudes 



We collect here the explicit values of the non-universal amplitudes appearing in (|3.ip - (|3.2p . 
which correspond to the normalized one particle/hole form factors (I3.5p - ()3.9p . They can be 
easily obtained from the expressions ()A.12p - ()A.2ip and read 



r q I 



2pl B[Fl^]e 



J[FlJ(q)-J[FlJ{-q) 



^ [2(zp'(g)] 



[Fl,{q) + l? + [Fl^{-q) + l? 

X 5? -^^-2 



jp sm 

2pl B[F-''] e^^''r]i-q)-J\FP\(q) 



T^{l + Fl^{q))T\l + Fl^{-q)) 



dz + TTFl^i-q) 



(5=0 



^ [2qp'{q)] 



{FlM-^]^ + \Fll-~q)~l? 



X dl siv? 



T^{l-F-'^{q))T\l-F-'^{-q)) 



dz + TTF-^iq) 



(A.28) 



/3=0 



in terms of the shift functions F-^^(A) = 4>{X, ^q) — (f'i^i ^q) = =tZ(A) =F Ij 



■Ao|2 



^q 



p(Ao 



epF 



Ao - eq 



Ao + Q 



27T^p'{Xo)p'iq) 



X djsm^ \i^Z + TTF^g'>{€q) 



(A.29) 



(5=0 



in terms of the shift functions F^I* = (/)(A, eq) — 4'{X, Aq), with e = ±1; 



KaoI 



p(Ao) 



Ao - eg 



- Ao\-^^".('") [2,y(,)]-[^".('^)l^-[^".(^'^)l^-^^^^o(^^) 



g + Ao 

X 



27r2p'(Ao)p'(g) 



X 5fl sin^ 



J[FZ]{eq)-J[FZ]{Xo) -p2 



r^(l + eFl^(eg)) 



i^Z(Ao) + vrF,;«(Ao) 



(A.30) 



15=0 



in terms of the shift functions F^^^ = (j){X,XQ) — (j){X,eq), with e = ±1. We recall that the 

functionals B and J are respectively given by (|X22]) and ([Al25l) . and that the expressions of 
the dressing functions are obtained through (|XT2l) - ([XJ3]) . 

Remark A. 4. For generic parameters, the /3-derivatives will apply directly to the factor (1 — k)^ 
of Qi (see remark lA.Sp . However, in the free fermion point c = +oo, they will apply on the sinus 
squared, since in that case the corresponding shift functions vanish (and Qi = 1). There may 
exist other particular values of the parameters for which the shift function becomes an integer 
but the overall expressions are anyway smooth (we may then have to apply the /3-derivatives to 
other factors such as the Barnes functions). 
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B Study of a generalized free fermionic generating function 

In this appendix, we consider the case of a generahzed free-fermionic model for which the 
effective generating function (j5.3p can be represented as a finite-size determinant. Namely, we 
suppose here that the shift function does not depend on the position of the roots parametrizing 
the corresponding excited state and that the dressing function Qn is separated (which happens 
in particular at the free fermion point c = +00 of the NLS model). 

More precisely, for a given counting function ^ (for example the ground state counting 
function of the NLS model), which defines a set of real parameters Aj, j = 1, . . . , by ^(Aj) = 
j/L, we consider the functional 



N 



N 

n 



sin^ [7rz^(Aj^ 



h<-<eN 



fi E^f^i,) fi TT^L^ e'(A,) eAj^e,) 



det 



1 



N 



where E__^ and u are holomorphic functions on some open neighborhood of 
expression should be understood as follows: 



(B.l) 



The above 



• the function S^^, is defined by = ^ — L u; 

• the multiple sum runs through all the possible choices of A^-tuples of integers £1 < ■ ■ ■ < £]\f 
belonging to the set ;Sl = {j G Z | —wl < j < wl}, where wl ~ L^'^'^, e > 0, is some 
cut-off which goes to +00 with L. 

• for a given set of integers Ij, the parameters are obtained as the pre-image of ij/L by 
the function This means, in particular, that they depend on v. 

Moreover, we restrict our study those functions that make (^j, a "good" counting function (i.e. 
ensuring a one-to-one correspondence between the rapidities and the integers). 

Proposition B.l. The functional X^v [i', E'^'j can be recast as the following finite size determi- 
nant: 



Xn [v, E'i] = det 



5jk + 



^e'(Afc) 

The corresponding finite-size kernel is given as 



(B.2) 



(B.3) 



where 



Ef\x)=iE^{X) { 



Bl 



— — — y + cot[7ri/(A)] + II [v, E_ \ (A) 



In this expression, the last term I/, [z^, i^^^] (A) corresponds to the integral 



lL[i^,EZ^]iX) 



dz EZ^{z) 



1 



27r z — X I 



+ 



dz Ez\z) 



1 



27r Z — X Q2iTTL(,^{z) _ I 



(B.4) 



(B.5) 
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Finally, the integration endpoints Al and Bl are such that L(^^{Al) = —wl—1/2 and L^^{Bl) = 
wl + 1/2, and '^1-/4. are some oriented contours, included in the joint domain of holomorphy 
of v and EZ^ ■ Theses contours join the points Al and Bl through the upper/lower half plane 
respectively, it has been depicted on Fig. [21 



Al 



Bl 



Figure 1: Contour "^-t 



t/i- 



Proof — The summand in (jB.ip is a symmetric function of the N summation variables ji^. 
which vanishes on the diagonals ij = ik, j ^ k. Therefore, we can replace the summation over 
the fundamental simplex £1 < • • • < ^tv in the A'^*'^ power cartesian product by a summation 
over the whole space Bf^ , provided that we divide the result by A^!. The summation domain 
being now symmetric, we can invoke the antisymmetry of the determinants so as to replace 
one of the Cauchy determinants by N\ times the product of its diagonal entries. This last 
operation produces a separation of variables, which enables us to recast the result into a single 
determinant by introducing the sum over £j into the j^^ line of the determinant: 



XM[iy,E^_ 



N 

n 



with 



jk 



4sin^[7ri/(Aj)] 
El{Xj).Ez'{ne) 



detvM, 



(B.6) 



EliX,) . Ez\f.e) 



+ (i-5,fc)y 



/i£ — Xj fj,£ — Xk 



(B.7) 



The above discrete sums can be expressed in terms of Hilbert transforms (plus corrections that 
vanish in the L — t- +00 limit). Indeed, 



EzHl^i 



E- - 



-iEZ^{X) 
-iEZ^{X) 

g2i7rLg^(A) _ Y 



+ 



EZ^{z)dz 



I 27r(z - A)(e2»'^^€-(^) - 1) 



EZ^{z)dz 
2tt{z - X) 



+ Il[v,EZ^]{X) 



EZ\z)dz EZ\X) 
2tt{z - A) 2 



cot[KLUX)]+lL[iy,EZ^]{X). (B.8) 
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Differentiating the above expression with respect to A leads to 



2sm'^ [7rL^i,(A)] 

By applying (fB^Sjl . dHU) in ([RT]) . and using the fact that ^u{Xj) = L^{Xj) -i^{Xj) = j-i^{Xj), 
we obtain 

M =5 g^(A,) j^-(A,) j^f (A,) ^-(A,) - (Afc) E-(A,-) 

4sin2[^i/(A,)] ^-(Afc) 2mL(Aj - A^) 

_ ^-(A,) ^^(A;.) f , I ^^^HA,,A,) 1 

i?_(Afc) 4sin[7rz.(A,)]sin[7rKAfe)]\ ^''^ ^^'(Afe) J' ^' ^ 

in which is given by (IB.4p . Note that we recovered the derivative of the A- type counting 
function ^ thanks to the identity = ^ — L~^v. This ends the proof of the proposition. ■ 

Let us now suppose that, in the thermodynamic limit L,N ^ +00 with N/L D [D 
finite), the roots Xj condensate on some symmetric interval [—q,q\ of the real axis with the 
density p[X) = lim^^jv-i-oo C'(A) of the NLS model. Let us moreover suppose that, in this limit, 
the counting function takes the form ()2.17p in terms of the dressed momentum p{X) of the 
NLS model. We demand in addition that log EZ^ has at most a polynomial growth in A when 
3fi(A) — )• ±00. Then the thermodynamic limit of ^^^[z^, -E^] is well defined and given by the 
following Fredholm determinant: 

lim XTvk^^l =detr_„„i[/ + y] (B.ll) 
with kernel 

Here, the function i?+(A) reads 

E,iX) = ^E_iX) I / + ^cot[vr.(A)]| . (B.13) 

Indeed, it follows from the form (|2.17p of the counting function that Al and Bl tend 
respectively to —00 and +00 in the thermodynamic limit. Moreover, it can easily be shown that 
Il\u,E'^\ = 0(L^^). The thermodynamic limit of Xi^[u,E'^\ is therefore a direct consequence 
of the fact that detAr[(5fc^ + o(L~"'^)] — )• 1, whenever the o(L~^) symbol is uniform in the entries. 



C Functional Translation operator 

In this appendix, we provide a heuristic approach to the notion of functional translation operator 
that we use in Section to express our highly coupled form factor series in terms of a decoupled 
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one. As we will see, such an object is in fact a convenient tool to manipulate generalized multi- 
dimensional Lagrange series (see Appendix C of |l5l)- We refer to [59j for a more explicit and 
rigorous construction. 

The notion of one-dimensional translation operator, which acts on some liolomorphic func- 
tion g as 



■g{uj) 



tj=0 



rv 



n>0 



(C.l) 



can easily be extended to the multidimensional case: 

s 

By analogy, we can thus define a translation operator acting on functionals as 

6 



(C.2) 



doo 7(0;) 



U[r] 



T = 



(C.3) 



Here, the integral is taken on a contour C (typically, an interval of the real axis), and f7 is a 
functional acting on holomorphic functions defined on a neighborhood U oi C. 

Formula (|C.3p can be understood as the result of the action of finite dimensional transla- 
tion operators on a finite dimensional approximation of U. More precisely, let us consider a 
discretization ti, . . . , of C and a collection of holomorphic functions Us{{zi}f^^) such that, for 
any holomorphic function 7, C/s({7(tj)}?i ) — > U[y]. Then, we can represent the translation 

s— >-+oo 

operator in terms of limits of the finite variable case: 



lim rrexp<^7(tp) 

p=l 



d 



Us{{Tm\ 



T=0 



(C.4) 



Such discretizations allow us to compute the action of more complicated translation opera- 
tors. Let us define, as in Section 15. 2| the operator ordered version : O : of some expression O 
containing functional derivatives of the type 5/5t{\) (such as in ()C.3P ) as being the expression 
where all functional derivative operators are placed on the left (in each term of the Taylor 
series expansion of O). Then, let us consider, for some functionals F and U, the following 
generalization of (10.3^ : 



/: = : expM F[r](A.)^3^d^J>[/[r] 
c 



+00 -If '1 



n=0 



6 



T=0 



.Sri,, [UnrK..)UM 



(C.5) 
(C.6) 



r=0 



After discretization one gets 



C = lim £, 



(C.7) 
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with 



p=i 



Tfe=0 



+00 s gm ( s 



- E n 

ni,...,ns £=i 
=0 



nil 



Y{r:^{{n}){r,)Us{{n}) 



p=i 



(C.8) 
(C.9) 



rk=0 



in which we have set = T{tk)- The last series is a multi-dimensional Lagrange series which 
can be computed as 



UsiMtk)}) 



dets [6jk - dr^Ts{{Ti}){t 



(C.IO) 



where i^itk), k = 1, . . . , s, are obtained as the solutions of the system z^(ifc) = Ts{{i'(t£)}){tk) 
k = 1, . . . , s. The s — )• +00 limit can be taken, at least formally. It yields 



detc 



:r[r](A) 



(C.ll) 



<5r(/i) 

in which v is the solution to the equation 

K/^) = rM(^). 



(C.12) 



D The Master equation issued-like series representation 



In this appendix we obtain, starting from (j5.12p - (|5.14p . an alternative series representation for 
Q'}f{x,t)es- This representation is in the spirit of the series obtained in [45j for the equal-time 
correlation functions. 



D.l The finite-size series 

Expanding the determinant in (|5.12p into its finite Fredholm series (this expansion is an imme- 
diate consequence of the Laplace expansion for determinants), we obtain: 



j=l ^-(^j) n=0 h,...,in=ls=l i^is) 



r=0' 
uj=0 

(D.l) 



m which f(^)(A,/i) is given by (I5.13p . It is convenient to use the analytic properties of the 
operator valued mappings E^_^^ and so as to re-cast this finite-size kernel V^^^ into a more 
compact contour integral representation: 



i>(-^) (A, ^) = 4 sin[7rz^^(A)] sin[7rz^^ (/i)] E_ (A) E_ (fi) 



dz E^_^\z)Ez\z) 
{2mY {z-\){z-^i) ' 



(D.2) 
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where surrounds [— (see Fig. [2]), A,/i G [— and 



2tt 



y 



EZ\y). 



(D.3) 



In this last expression the contour ^^-^^ = 'loj^ VJ'^qVJ'io^VJ'io^ consists of a union of four contours. 
Two of them, "^-^/i, are as depicted in Fig. [H and the remaining two, 'rfj^^U^g, appear in Fig. O 
In particular, the contour 'ifq encircles the contour 'Wq of ()D.2p . which enables us to recast the 
integral kernel ^(^^(A, jj,) in a more compact form. The integrand in ()D.3P involves the function 



(y) 



+ 



- + — 
1 1 



1^4 (y) 



(D.4) 



(L) 



Al 



— ^ 

r 



Bl 



Figure 2: Contours '^q and '^q. The endpoints Ai^ and Bi^ are such that Li^y^{^Ai) 

-WL - 1/2 and L^^^{Bl) = wl + 1/2. 



Factorizing the integrals over z ()D.2p out of the determinant in (jP.ip and then using the 
symmetry of the summand in order to reconstruct a second Cauchy determinant, we get 



N 



d"z 



'2i7r 



,2n 



-2i7r)' 



n 

s=l 



1 4 sin^ 


7rz^r(AiJ] 






1 • det^ 


1 




_ Zj — Xif, 



r=0 
w=0 



Since 



^2 _ ^-ixu-gr-du 



(D.5) 
(D.6) 



it remains to compute the effect of the two types of functional translation operators occuring 
in (ID.Sp . We start by taking the w-translation into account, namely 



a=l a=l 







uj=0 





{A,}f \{A.JLiU{2/s}:=i 



N 



(D.7) 



Above, we have first computed the action of the translation operators according to (15. 4p . then 
used the equality (|4.13p . We should now compute the action of the second set of translation 
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operators involving Qt- Since the parameters fii = /1£[i't] = t= 1, . . . , A^, at which the 

are evaluated, are themselves functionals of r, we should take care of the operator ordering. 
The expression we have to compute is of the type: 



n— u Ji,...,j„— i ^f^-^ Ij, 

,w ({ysYl 



{^4 \{zs}'l 



, (D.8) 

r=0 



in which } is a smooth functional of I'j- that, moreover, depends on the integration variables 
{ys} and {zg}, and the functional r|j^}[r](A) driving the functional translation reads 

N n 

riijT]{X) = Y,m^,fii['^r])-Hx,Xe)] + Y.[4>{x,KJ-4>{\ys)]. (d.9) 

e=i s=i 

According to Appendix O (see (jC.Sp ). computing the action of a functional translation operator 
with weigths r|j^}[r](A) as above amounts to summing up a multi-dimensional Lagrange series. 

(n) 

In our case, the result is obtained by replacing the functional argument i^r of TZ}.\ by the 
function z/j, in which r is the solution to the non-linear functional equation 'c(t) = rjj^j.[r](t), and 
then by dividing the obtained expression by the corresponding functional Jacobian J{j^ j[r^] = 
detM [/ - ,5%}[r](A)/<5r(A)] |^^^ (see ([0311)). 

Finally, Q'j^{x,t)eS admits the following representation: 



ra— u ii,...,i„— i ^{L) 



X 



.=1 yLi'{Xi^){ys- Zs) 



X det^ 



n 



In (ID.lOp . we have set 



N 



u\'^y{Xa]AUa[^\]Ays]) = Y.^u{X,)-u{H^])] + Y.[u{y,)-u{Xa (D.ll) 

and explicitly insisted on the fact that Jin = ^ = 1, . . . , A^, are now functionals of the 

solution r' to r'(t) = r{j^}[t'](t). 

D.2 Taking the thermodynamic hmit 

For finite A^, (jP.lOp gives a rather implicit representation. We are however interested in com- 
puting the thermodynamic limit L,N ^ +oo of this expression (in fact it is only in this limit 
that the effective series (ID.lOp is supposed to coincide with the original form factor series (14. BP ). 
In this limit of interest (|D.10p can be considerably simplified. 
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In particular, the non-linear functional equation 'i{t) = T ^^^y[t]{t) for r turns into a linear 



integral equation for r(t), with r(t) — > r(t): 



x{t) = v^{t) — i 



2n 



d^Ht, A) z.,(A) dA + A,J - cPit, vs)]. 



(D.12) 



s=l 



Indeed, in the thermodynamic limit, the Bethe roots Xj for the ground state condensate on 
[— g, q] with the density /9(A), whereas the parameters (defined in terms of the counting 

function are shifted with respect to the Aj's according to: 



Lp{Xj) 



j = l,...,N. 



(D.13) 



Above, we have denoted the thermodynamic limit of the function z/j by f^. The linear integral 
equation (jD.12p is readily solved as soon as one observes that the derivative f?2</> of the dressed 
phase (I2.12P with respect to its second variable is actually related to the resolvent R of the Lieb 
kernel {d2^{t, A) = -R{t, A), with (/ - K/2tt){I + R) = I): 



a 1 " 



(D.14) 



It also follows from ()D.12p that the Jacobian of the transformation (see (jClip ) is simply given 

by 



J = det[_g^g] [/ - d2(t)] = detrj^ 



27r. 



(D.15) 



The expression (ID.14P for the thermodynamic limit of i^{t) along with the estimations for 
the shift (jP.lSP of the parameters fij [z^] with respect to the Aj's allows one to compute the 

thermodynamic limit of Z^j^ j({Aa}, {/ia}, Namely, 



Uf^^{{Xa},{-^a},{ys}) 



N,L—i-oo 



- u'{x)v,{x)dX + Y,<ya)-u{K^) 

-q ^^=1 

= -i^ + E[no(ya) - no(A,J]. (D.16) 



a=l 



To obtain this limit, we have used the definitions (I2.15p -( f2.16p of e and p so as to re-cast the 
expressions involving integrals of u = p — te/x in terms of uq = po — t£o/x. 

The thermodynamic limit of can be computed along the lines of |43] . Indeed, one has 



F X 



{A,}f U{yj:=i 



z^r(A) . 



The latter implies that 

N 

\ 

io — Xj + iec N,L^oo 



(D.17) 



■j-j- id - flj + iec _ ^-2inClF]{uj+iec) ^ ^ 



for eG{±l,0} 



(D.18) 
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By applying identity (jP.lSp to the various products entering in the definition of ^^r, one can 
convince oneself that 



det[_,,] [/ - K/2vr] n4sin^[vr..(AJ] ^.v (^^^'^^ ^ ^^^^!.^^^^^^ 



s=l 



N.L^OD 



n 



{VaYl 



{ys}1 



(D.19) 



Above, we have set 



{^i},=i 



N 

n 



^Ji =F (// - Za) 

and, agreeing that the auxiliary arguments of V± are undercurrent by those of 



(D.20) 



J' n 



(1 - k)' 



{As}? 

{ys}^J det[/-ii72vr] 11 



n 



jk J 



n 



iVa - Afo - ic) {yb - Xg- ic) 
^ iVa - yb- ic) {Xb - Xa- ic) 



, (D.21) 



with 



V,. 



(8) 



( k y^) j^^^^^ _ ^^-^ _ j^^f^Q _ ^^-^ 



kl 



n (Afc - A. 

s=l 

^k 



vi'iXk)-KVr'{Xk) 



(D.22) 



It only remains to deal with the y-integrations over In this limit, the integrals over 

"^t/i S^^^ vanishing contributions and hence the y-type integration contour boils down to the 
contour ^ = "^^""^ U ^ with a weight function that now reads 



f{y,'^x (y)) = Voc) (y) 



Q-2i-KU,:{y) _ -y 

to the L — )• + 
Therefore, in the thermodynamic limit, Q^(x,t)efr 



(D.23) 



The contour corresponds to the L — )• +oo limit of the contour "^^^^ depicted on Fig. [21 



Q'^{x,t), with 



n«r -V" (~^)" / if d"z /- 

^'''^ " t^o (^J (^i i^^-r t\ 



d"y 



gix[Mo(2ya)-Mo(Aa)] 



n 

s=l 



f{ys,Mys)) 



1 - K- 



{Aa} 



detf 



Afc 



{Aa} 

{ya} 



(D.24) 



Here, we remind that i/^ is given by (ID.14P and is a function of {ys} and {Aj}. 
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D.3 The equal-time case 

We now show that in the equal-time case, one recovers the series obtained from the master 
equation in |45) . This stems from the fact that, at t = 0, the y-integrals over 'j^^^ do not 
contribute. Indeed, when t = 0, we can deform the contour '^^^ into M + i6, with 6 > and 
smah but such that the hne M + i(5 is above '^q. 

In order to prove this assertion, we first build on the symmetries of the integrand in (|D.24p 
so as to split the y- integrations into those along M + z5 (i/a, a = 1, . . . ,k) and those along 'Wq (ija, 
a = k + 1, . . . , n), with k = 1, . . . , n. The integrals over '^q can then be computed by residues. 
One eventually obtains 



(-ir 



n=0 k=0 



d"A 



nlkl{n-k)l J {2mY J {'2iTr)'^ J {2mY 



d"z 



X deiz 



n 

a=k+l 
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n 
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s=l 


V ( 


{ya}\yj{Za}UJ_ 


Zj — Afc 


n 

n 

a=k+l 




Za 





^ n, 



n 

a=l 



Jx[ya-Xa] 

Va - Za 



{ya}\^{Za}l+^ 



(D.25) 



When considered as a function of yi , . . . , , the integrand in (ID.25P has no poles (or even 
other singularities) in the half-planes 9 (y^) > 5. Indeed, no poles can arize from the 9j 
dependent terms in in as much as this function does not depend on Oj, cf |45| . Also the 
potential singularities of the determinants at the zeroes of 1 — K(y_/V+)(Aa) are only apparent 
due to the presence of the pre-factors — K~^(y+/I^_)(As)][l — K(y_/y+)(As)] distributed in 
between (|D.2ip and (|D.25p . The potential poles at y^ = A^ it ic introduced by these pre-factors 
are cancelled by the zeroes of the double product in the last line of (ID.2ip . In its turn, the 
poles in the upper-half plane at ya = Zg + ic, with a = 1, . . . ,k and s = k + 1, . . . ,n, that are 
introduced by the double product in ()D.2ip . are compensated^ by the same poles present in 
— K,{V-/V-j-){zs)]. Therefore the only singularities in the y variables correspond to the 

zeroes of 



{Aa}^ 
{ya}1^{Za}l+^ 



(D.26) 



However, one can always squeeze the z-integration contours in (ID.25P so that ^(za) 
a = 1, . . . , n. In such a situation, it is readily seen that for y^ G M -|- i5, 6 > 0, one has 



1^^^ A., 



> 1 



(D.27) 



In other words, the aforementioned equation has no solutions in the upper half-plane. This lack 
of singularities allows one to simultaneously deform the y-integration contours from M + i(5 to 
M + iA, with A > and as large as desired. Due to the presence of the oscillatory factors e*^^", 
these integrals will contribute as e~^^. Therefore, by sending A — )• +oo, we see that these 
contributions vanish. 



^ The poles at ya 
for our purposes 



Zs — ic lying in the lower half-plane are not explicitly compensated but this is irrelevant 
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We have thus proven that, in the t 
for Q'^{x,0) boils down to 



case, the series of multiple integral representation 



Q-(^,0) = y^ /-^ /^!^fre-(^-M det 
^ ' ^ ^ (n!)2 J (2i7r)" J (2i7r)" JJ- 



{Aa} 



n 







{Aa}A 









. (D.28) 



In order to prove our assertion, it remains to show that the last line in (ID.28P does not contribute 
once that the 2;-integrals are taken. In virtue of the symmetry of the integrand, one can 
replace one of the Cauchy determinants by n! times the products of its diagonal entries. By 
expanding the remaining Cauchy determinant into a sum over permutations, one sees that 
one has to compute double poles at Zp^ = Ap^. with pi < • • • < pi with pi G {1, . . . n} and 
1 < i < n. All other poles will be simple, leading to the equality between the two sets 
{za}aytpf, = {Aa}a7^pif These double poles will produce first order derivatives with respect to 
^Pk -Zp,. = Ap^.. Therefore, the effect of these double poles can be taken into account by 
setting^ Za = Xa for a ^ pk, k = 1, . . . and = Ap^, + ep^, and then taking the first order 
epj, -derivatives of the obtained expression at gp^. =0. As a consequence, only the linear in each 
epj. order of the integrand will contribute. However, under such a substitution, one can readily 
convince oneself that 



n 1 ^ 

n— 

«=1 1 - K 



V- 




{Aa} 




(zs 


{Aa} 



IK 



4^ 



The linear order in the ep^. 's vanishes. Therefore, we are led to 



Zj — Afc 



J-„,(|^"jj. (D.29) 



Once upon taking the complex conjugate we recover, word for word, the series obtained in |45| . 
Remark D.l. Note that our conventions correspond to x i— t- —x with respect to the work 



E Controlling sub-leading corrections: the Natte series repre- 
sentation 

The Natte series representation for the Fredholm determinant det [/ + V] is obtained |6Q) from 
a specific representation of the solution to the Riemann-Hilbert problem associated with the 
integrable integral operator I + V . As this particular representation for the solution of the 
Riemann-Hilbert problem is obtained by a series of contour deformations on the so-called initial 
Riemann-Hilbert problem associated with the integral operator I + V , it can be seen that the 

i^The integrand is symmetric with respect to the integration variables Za and A^. 
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Natte series stems from a certain number of algebraic transformation^^ carried out on the 
initial Fredholm series for the determinant. This fact allows one, at least on the formal ground, 
to use this Natte series representation in (|6.5p . The rigorous justification of the possibility to 
use the Natte series is given in |59| . 



E.l The formula for the remainder 

In fact, the Natte series corresponds to a representation of the remainder TZxl^, u, g]. The latter 
is expressed as a series of multiple integrals. The integrands appearing in this series have good 
properties with respect to the large-x limit. 

More precisely, the Natte series for the sine kernel (16. 2p takes the following form: 



The second sum appearing above runs through all the n-tuples k belonging to 

(E.2) 

For each element k of /C„, one defines the associated set of triplets J(k) (with cardinal n), 

J(k) = \t={ti,t2,t3) , ti G [1; nl, t2 G [1; /ctj, t3 G [1; till , (E.3) 



/C„ = < k = (/ci , . . . , kn) : ka € ^ ,a = 1, . . . ,n and such that skg = n> 

^ s=l J 



which provides a convenient way of labelling sets of n variables. The third sum runs through 
all the possible choices of elements belonging to the set 

£:„(k) = |{et}t6j(k) : ete{±l,0} and X] et = V(ti,t2) G fl ; x [1 ; . 

In other words, i?n(k) consists of sets of n parameters et E {±1,0}, indexed by triplets t = 
(^15^2,^3) G ^(k) and subject to summation constraints. Finally, there is an n-fold integral 
appearing in the n}^ summand of (lE.ip . The integration variables zt are, again, indexed by the 
triplets t = (ti, ^2; ^3) of J(k). The contours of integration '^{ej} depend on the set {et} G <?n(k). 
They are realized as n-fold Cartesian products of one-dimensional compact curves corresponding 
to various deformations of M. 

The integrand Hn-xH^t} ] {zt})['^] is a smooth functional of which is also a function of 
the integration variables Z{. This functional depends on the choice of the parameters {et} from 
iSn(k) and on x. It is a holomorphic function of the integration variables {zi} belonging to some 
open neighborhood of the integration contour "^{et}- 

The Natte series converges for x large enough in as much as, for n large. 



/An _ 
^/7„,({et},{.t})M n 



<e.(Si)"". (EA) 



X 



^Such as contour deformations or algebraic summations. 
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where ci and C2 are some constants depending on the values taken by u and g in some small 
neighborhood of M and on the values taken by on a small neighborhood of [—q^q\- The 
constant C3 > only depends on v. 

Finally, one has Hi-^x = 0{x~°°) and, for n > 2, 



i^„;x.({6t},W)M=0(x— ^ f 

6=0 p=0^=b_[|] 



X2[u{-q)\ \ \ ^2[u{q)+u(-q)\ 



x-i^. i7-f''({6t},{zt})M, (E.5) 



where, for all admissible values of m,b,p, 

^nT'({^t} ' = 0(x"--^-^) (E.6) 

uniformly on the integration contour. Above, we have set 

w = 2sup!^\R[u{z) -u{Tq)]\ : \z - Tq\ < 6 , t = ±^ , (E.7) 

where 6 > can be taken as small as desired. These functions H^^^'i'^ ({^t} ;{-2t}) [^] are such 
that their asymptotic expansion into inverse powers of x possesses poles that are encircled by the 
part of the contour ^{e^} producing algebraic (in x large enough) contribution to the integral. 
As a consequence, the n-fold integration occurring in the n*^ term of the series can produce z 
derivatives. It was shown in |60| that the total order of these derivative is at most n. Once that 
these derivatives are taken, the estimates in ()E.6P change from 0{x'^'^~^^'^) to O^x"^^'^ log" x). 

E.2 The Natte series for the determinant 

We now use the Natte series representation (jE.ip for the remainder TZx\i',u,'g\ to prove that 
the subleading corrections appearing in the large-x asymptotic behavior of the (non operator- 
dependent) Fredholm determinant and included in TZxli', u,g] still remain corrections with re- 
spect to the leading terms once that one computes the action of the functional translation 
operators. By expanding the Fredholm determinant in (16. 5p into its Natte series, we obtain 



[ e=±l 

+ \Y1 e*"^["(^»)-"(^«)l+"[^(^o)-^(^«)l ^^''"^[z^^,^] Bx[iyr,u] 



e=±l 



n>l K.„ £„{k)cJ , ^ teJ(k) J 



w=0 

T = 



(E.8) 



Note that in (jE.Sp . we have already simplified the two exponents. Also, we remind that a = 1 
in the space-like regime and a = — 1 in the time-like regime. The action of the translation 
operators occuring in the first two lines has already been computed in Section 16.31 In order to 
compute this action in the last line of (lE.Sp . it is convenient to introduce the notations 

a{^^l = #{t : et = 1} and {zf} = {zi : t G J(k) such that ± et > 0} . (E.9) 
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It then readily follows that 



n[e--(-)]-^M| =S.i|^t|), (E.IO) 

tGJ(k) """" ^^^t 

and 

: Yl [e^'^.(^.)] .F„;,({et},{zt})K] e,K,n] : |^^^= /7„;,({6t}, B.AF{,^y,u]. (E.ll) 
teJ(k) 



Note that, in order to obtain (lE.lOp . we have used the fact that et G {±1,0} for all t G J(k) 
and that these parameters are subject to the condition X]teJ{k) = 0. As a consequence, 
= Also we have set 

i^{e.}(A;{^t}) = = ^.(A,{zt}) = '-^^ - J2 (E.12) 

teJ(k) 

This leads to the following representation for the generating function: 

d"zt 



(2i7r)" 

.+1 



X i^n;x({et},{^t})[F{.,}] B4F^,,},u] . (E.13) 

Here, Qasymix,t) denotes the leading asymptotic part that is given in (I6.15P or by its similar 
expression in the time- like regime (without the 0((lnx)/x) remainders since these are now 
explicitly taken into account in the rest of the above equation). 

Note that the n-fold integration occuring in the n^^ term of the series can produce z deriva- 
tives whose total order is at most n. It thus follows from the representation for the functional 
Hn-x (jE.SP and from the form of the estimates ()E.6P that, indeed, the remainder produces 
corrections of the form written in (I6.15P . 
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